Age and size compositions

Fishery age composition

Code

library(dplyr)
library(ggplot2)
library(plotly)
theme_set(theme_bw())
# data ----
source(here::here('r', 'bridge_data.r'))
data.frame(age = ages,
             year = readRDS(here::here('data', 'rtmb', 'report1.rds'))$fish_age_pred,
             groups = 'rtmb') -> df

names(df) <- c('age', fish_age_yrs, 'groups')
read.csv(here::here('data', 'admb', 'fac.csv')) %>% 
  mutate(groups = ifelse(groups=='pred', 'admb', groups)) %>% 
  bind_rows(df %>% 
             tidyr::pivot_longer(-c(age, groups), names_to = 'year') %>% 
      mutate(year = as.numeric(year))) %>% 
  ggplot(aes(age, value, color = groups, linetype=groups)) + 
  geom_point() + 
  geom_line() +
  facet_wrap(~year) +
  scico::scale_color_scico_d(palette = 'roma', end = 0.8) +
  scale_linetype_manual(values = c(0,1,1), guide = guide_none()) -> f1

ggplotly(f1)

Survey age composition

Code

# data ----
data.frame(age = ages,
             year = readRDS(here::here('data', 'rtmb', 'report1.rds'))$srv_age_pred,
             groups = 'rtmb') -> df

names(df) <- c('age', srv_age_yrs, 'groups')
read.csv(here::here('data', 'admb', 'sac.csv')) %>% 
  mutate(groups = ifelse(groups=='pred', 'admb', groups)) %>% 
  bind_rows(df %>% 
             tidyr::pivot_longer(-c(age, groups), names_to = 'year') %>% 
      mutate(year = as.numeric(year))) %>% 
  ggplot(aes(age, value, color = groups, linetype=groups)) + 
  geom_point() + 
  geom_line() +
  facet_wrap(~year) +
  scico::scale_color_scico_d(palette = 'roma', end = 0.8) +
  scale_linetype_manual(values = c(0,1,1), guide = guide_none()) -> f1

ggplotly(f1)

Fishery size composition

Code

# data ----
data.frame(length = length_bins,
             year = readRDS(here::here('data', 'rtmb', 'report1.rds'))$fish_size_pred,
             groups = 'rtmb') -> df

names(df) <- c('length', fish_size_yrs, 'groups')
read.csv(here::here('data', 'admb', 'fsc.csv')) %>% 
  mutate(groups = ifelse(groups=='pred', 'admb', groups)) %>% 
  bind_rows(df %>% 
             tidyr::pivot_longer(-c(length, groups), names_to = 'year') %>% 
      mutate(year = as.numeric(year))) %>% 
  ggplot(aes(length, value, color = groups, linetype=groups)) + 
  geom_point() + 
  geom_line() +
  facet_wrap(~year) +
  scico::scale_color_scico_d(palette = 'roma', end = 0.8) +
  scale_linetype_manual(values = c(0,1,1), guide = guide_none()) -> f1

ggplotly(f1)